Superconductivity in zigzag CuO chains 
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Superconductivity has recently been discovered in Pr2Ba4Cu70i5_i with a maximum T c of about 
15K. Since the CuO planes in this material are believed to be insulating, it has been proposed 
that the superconductivity occurs in the double (or zigzag) CuO chain layer. On phenomenological 
grounds we propose a theoretical interpretation of the experimental results in terms of a new phase 
for the zigzag chain, labelled by C1S3. This phase has a gap in the relative charge mode and 
a partial gap in the relative spin mode. It has gapless uniform charge and spin excitations and 
can have a divergent superconducting susceptibility, even for repulsive interactions. A microscopic 
model for the zigzag CuO chain is proposed, and on the basis of density matrix renormalization 
group (DMRG) and bosonization studies, we adduce evidence that supports our proposal. 



Introduction - The discovery of high-T c superconduc- 
tivity has raised the possibility of new mechanisms of su- 
perconductivity in strongly correlated electron systems, 
which arc radically different from the well established 
BCS-Eliashberg mechanism. However, the lack of the- 
oretical techniques for obtaining well controlled solu- 
tions of even the simplest models of strongly correlated 
fermions in dimensions greater than one has seriously 
limited the understanding of such mechanisms. 

An important exception is one dimensional systems, 
where powerful analytical and numerical techniques are 
available. A one dimensional quantum system cannot 
have superconducting long-range order even at zero tem- 
perature, as a consequence of the Mcrmin- Wagner theo- 
rem. Nevertheless, such a system can have a supercon- 
ducting susceptibility which diverges as T — > 0. There- 
fore, an array of coupled one dimensional systems can 
become a "bulk" superconductor with a reasonably high 
T c even in the limit where the inter-chain coupling is very 
weak. This mechanism of superconductivity has been 
studied extensively^. It serves as one of the only well- 
established proofs of principle for superconductivity in a 
model with purely repulsive interactions. 

Among the candidates for an experimental realiza- 
tion for this mechanism are the organic conductors and 
the ladder compound Sr^-zCa^C^^^. Recently, su- 
perconductivity was discovered in P^Ba^uyOis-a (Pr- 
247) in a certain region of oxygen reduction 8, with a 
maximum T c of about This material is isostruc- 

tural with Y 2 Ba4Cu 7 15 _ < 5 (YBCO-247). However, the 
two materials display a dramatic difference in their elec- 
tronic behavior: the copper-oxide planes in YBCO are 
conducting, and are believed to play a crucial role in 
the high-T c superconductivity in this material. In Pr- 
247 (as well as in the closely related material Pr-248) the 
copper-oxide planes remain antiferromagnetic with large 
moments even upon doping. To the best of our knowl- 
edge, the conductivity of a single crystal of Pr-247 has 
not been measured yet. In Pr-248, the conductivity is 
strongly anisortopic in the a — b plane (with o- a /&b of up 
to 1000, where a is the chain direction&i^) . This strongly 



suggests that the copper-oxide planes in these materials 
are insulating, and that most of the conductivity occurs 
in the metallic copper-oxide chain layers. Since the dou- 
ble (or zigzag) chains are much more structurally robust 
than the single chains, they are much less disordered and 
therefore are expected to be better conductors. 

Assuming that the electrical conductivity in Pr-247 (as 
in Pr-248) comes mostly from the double chains it follows 
that the superconductivity must originate in these chains 
as well. While this assumption remains to be tested, ex- 
isting NQR experiments^ which measure the site-rcsolvcd 
spin-lattice relaxation rate (1/Ti) have shown that only 
the double chain copper nuclear spins show any sharp 
feature (a cusp) in their 1 jT\ (T) curves at the supercon- 
ducting transition temperature, supporting the idea that 
the superconductivity is intimately related to the double 
chains. This possibility can have important implications 
in other materials which also share the same zigzag CuO 
chain structure, such as YBCO-247, YBCO-248 and the 
ladder compounds. 

The purpose of this paper is to study the possibility of 
superconductivity in zigzag CuO chains. This problem 
has been addressed in several previous papers^ ' 10 ^ 1 . In 
Ref. 0, the CuO zigzag chain was treated in the weak 
coupling limit and in the Hartree-Fock approximation. 
Ref. |10| st udied a zigzag Hubbard ladder using DMRG. 
In Ref. ll[ weak coupling fluctuation exchange (FLEX) 
theory was used to study superconductivity in the zigzag 
chains. 

In this paper, we start from the experimental data, 
and analyze the experimental constraints on the zigzag 
chain superconductivity scenario. Then, we present a 
microscopic model for a single zigzag chain, which con- 
tains in our view an important piece of physics that has 
been omitted in previous studies, namely the oxygen p y 
orbitals and the ferromagnetic Hund's rule coupling on 
the oxygen sites.— The model is studied using numeri- 
cal density matrix renormalization group (DMRG) cal- 
culations, as well as analytic renormalization group (RG) 
and strong coupling methods. The solution of the model 
shows features that are in agreement with experiment, 
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namely the absence (or extremely small magnitude) of 
the spin gap and a tendency towards superconductivity. 



I. POSSIBLE PHASES OF THE ZIGZAG CHAIN 

Our main assumption is that the zigzag CuO chains 
in Pr-247 are weakly coupled, so the basic properties of 
the system at T > T c can be understood in terms of 
the properties of decoupled chains. Let us review briefly 
a mechanism of superconductivity in weakly coupled Id 
systems. A Id system typically has several types of fluc- 
tuating order which coexist with each other (for exam- 
ple, superconducting and spin or charge density wave 
orders). Even though none of these can truly become 
long-range ordered in the isolated Id system at generic 
(incommensurate) filling, the susceptibility of the system 
to these types for order can become large at low tem- 
perature. Then, treating the inter-chain coupling at the 
mean-field level, the critical temperature is determined 
by the Stoner criterion: Jx{T c ) = 1, where J is the inter- 
chain coupling and x{T) is the susceptibility to the type 
of order considered. Assuming that the inter-chain cou- 
pling constants of the various types of order arc all of 
comparable size, the type of order that is most likely to 
be selected is the one which has the largest (i.e. most 
divergent) susceptibility at low temperature. 

A single-component Id electron system with repul- 
sive interactions can generically be described at low en- 
ergies as a Luttinger liquid with one gapless spin and 
one gapless charge mode. Such a system is usually a 
poor superconductor, with a superconducting suscepti- 
bility xsc ~ T x l Kc ~ x , where K c is the charge Luttinger 
parameter. Since typically K c < 1, the superconducting 
susceptibility is non-divergent while the CDW and SDW 
susceptibilities both diverge. 

Multi-component Id systems are more promising in 
this respect. The problem of two coupled Luttinger liq- 
uids has been studied extensively, and it is found to be in 
a "Luther- Emery" phase over a wide range of parameters. 
In this phase, only the total charge mode remains gap- 
less while all other spin and charge modes are gapped. 
It is therefore similar to the phase of a single compo- 
nent system with a spin gap caused by attractive inter- 
actions, although the identity of the physical correlation 
functions in the two systems is somewhat different. If 
we label the Id phases as C ra S m , where n and m are the 
numbers of gapless charge and spin modes respectively^, 
the Luther- Emery phase is CiSo- Since the entire spin 
sector is gapped, the dominant fluctuations in this sys- 
tem are superconducting fluctuations and 4k p CDW fluc- 
tuations, with susceptibilities \sc ~ T 1 ' 2K +' C ~ 2 and 
XCDW,4k F ~ y 2A: +,c-2 respectively. (Here, K +)C is the 
total charge mode Luttinger parameter.) Near half fill- 
ing, it has been showni^ that K +tC ~ 1 and therefore 
superconducting fluctuations always dominate. The in- 
terpretation of superconductivity in Pr-247 as deriving 
from a Luther- Emery phase was proposed in Ref . Q . 



An important observation is that in this scenario, the 
spin gap A s must be significantly larger than T c . At tem- 
peratures higher than A s , the spin modes are essentially 
gapless, and the exponents controlling the susceptibilities 
would cross over to new exponents of a gapless phase. In 
the regime T c < T < A s "pseudogap" behavior should 
be observed: a gap appears in the spectrum, but this gap 
is not associated with any type of long-range order. 

However, NQR measurements of the spin-lattice relax- 
ation rate \/T\ above T c show no evidence of gapped (i.e. 
thermally activated) behavior^. Assuming that the re- 
laxation process is mostly due to local magnetic field fluc- 
tuations (rather than electric field gradient fluctuations) , 
this implies that there is no spin gap in the zigzag chains. 
The 1/Ti signal seems to follow a non-trivial power law 
as a function of temperature, which implies a power law 
decay of spin-spin correlations. For a discussion of the 
interpretation of 1/Ti NMR measurements, see [27| . 

The preceding analysis leads us to the phenomenolog- 
ical idea that the zigzag chain is in a phase with gapless 
spin excitations. Two uncoupled chains are in a C2S2 
phase. Inter-chain interactions can gap some of these 
modes. For example, we can consider a C1S2 or a C1S1 
phase. In these cases, the spin-spin correlations behave as 
a power law (since the total spin mode is gapless) . How- 
ever, some of the relative spin and charge modes (the 
modes that correspond to fluctuations of the charge of 
the two chains relative to each other) are gapped, which 
can enhance superconducting (as well as other) suscep- 
tibilities. In the rest of this paper, we will advocate an 
explanation of superconductivity in Pr-247 in terms of 
such a phase, whose existence will be established in a 
microscopic model for the zigzag chain. 

The most likely candidate phase that emerges from our 
analysis is what we will call a C1S3 phase, in which the 
relative charge mode and "half of the relative spin mode 
are gapped. (We give a more precise definition for these 
terms in Appendix [All The resulting susceptibilities in 
this phase are 

Xsc ~ TVa* + ..-f (1) 
XcDWAk F ~ T 2K +>°~ 2 (2) 

XCDW,2k F ~ XSDW,2k F ~ T K+ '"^ 2 ~^ (3) 

Note that the superconducting susceptibility is divergent 
for K + c > 0.4, and is the dominant one when K +c > 1. 

We comment that the proposed C1S3 phase is likely 
not to be a true T = phase: it will be ultimately un- 
stable to formation of a C1S0 phase. Nevertheless, if the 
gap in the total spin sector is small compared to T c , the 
physics above T c is best described in terms of a CiS|. 
Such a hierarchy of gaps is found to emerges in the mi- 
croscopic model that will be studied here (see Section 

Finally, the Luttinger exponent K +tC that controls the 
low-energy properties of the system depends on micro- 
scopic details, and is therefore difficult to estimate theo- 
retically. However, since the power with which the spin 
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FIG. 1: The geometry of the zigzag CuO chain. Squares mark 
copper atoms and circles are oxygen atoms. Also shown are a 
copper d x 2_ y 2 orbitals and an oxygen p x orbitals. The dotted 
box marks the unit cell. 



susceptibility depends on temperature also depends on 



-K+,o it can be extracted from a measurement of this 
susceptibility. From the NQR measurement of \jT\T ~ 
Xsdw we extract the value K +tC w 3/2. Note that this 
corresponds to a regime of dominant superconducting 
fluctuations with xsc ~ T~ n / 12 . 



II. THE MODEL 

We now present a model for a single zigzag CuO chain. 
The geometry of the zigzag chain is shown in Fig. [TJ Of 
this structure, we will ignore the outer oxygens, leaving 
two copper and two oxygen atoms per unit cell. The 
relevant orbitals are the 3d x 2_ y 2 on the Cu sites, and the 
2p x and 2p y orbitals on the O sites. We assume that 
the on-site Coulomb repulsion is large on both Cu and 
O sites, so that doubly occupied states can be projected 
out, leaving an effective magnetic exchange interaction. 
The Hamiltonian for holes can then be written as 



H = H o + H e 



(4) 



Here Hq is the hopping Hamiltonian and H ex is the mag- 
netic spin exchange Hamiltonian. Let us first consider the 
terms in Hq. A 3d Cu orbital can hybridize strongly with 
the p x orbital of its neighbor O in the same chain, or a p v 
orbital of the neighboring oxygen in the opposite chain. 
However, it cannot hybridize with the nearest p x orbital 
in the opposite chain or a p y orbital in the same chain 
due to symmetry. Further neighbor hopping matrix ele- 
ments (such as a direct oxygen-oxygen hopping t pp shown 
in Fig. [T|) will be neglected, except in section HV Bl where 
their effect on the main results will be examined. For 
simplicity, we also also project out the p y orbitals, since 
they essentially "belong" to the neighbor Cu d orbital, 
in the sense that Hq only connects a p y to that d orbital. 
The effect of the p y orbitals will be reintroduced as an 
effective interaction in H ex . Taking these considerations 



into account, Hq is simply 

Hq = t pd J2 P (4<7,J»s,i+l,<7," - Pl,i,<T,oc d i,<r,cx + h - C ) P 

(5) 



Here d ij(T , Q and p x . 



are hole annihilation operators 



in 3d Cu and 2p x O orbitals respectively, P is a pro- 
jection operator that imposes the no-double occupancy 
constraint, i is the unit cell index, a =T,j is the spin 
index and a = 1, 2 is the chain index. The oxygen sites 
have an on-site potential e > 0. 



Next, we consider H e 



H„ 



This has the form 
1 



; n d. : 



,a.Tld,i-\-l,Oi 



(6) 



Here S d = J2a,a' 4^'4' where a = (a x , <r y , a z ) are 
Pauli matrices, rid = and similar definitions for 

S Px , n Px . J\ > is the usual antiferromagnetic superex- 
change interaction. J2, however, has a qualitatively dif- 
ferent physical origin: it comes from projecting out the 
state with a doubly occupied oxygen site with one p x 
state and one p y state occupied by holes. Since the two 
holes belong to different orbitals, they are subject to a 
ferromagnetic Hund's rule interaction^. Therefore the 
effective exchange interaction is ferromagnetic (J2 > 0) 
in this case. (Note the minus sign in Eq. [SI) 

The model we have derived for the zigzag chain has 
several new features that distinguish it from the well- 
studied Hubbard ladder model. Firstly, as we have seen, 
the inter-chain hopping is small relative to the intra-chain 
bandwidth. (It arises only from further neighbor hop- 
pings, such as the 0-0 hopping.) The coupling between 
the two chains thus comes mostly from the electron- 
electron interactions^. Secondly, the inter-chain effec- 
tive exchange interaction is ferromagnetic. 



III. HALF FILLING 

It is instructive to consider the case of half filling in the 
double chain, which we consider as the parent insulating 
state. In this case (i.e. one hole per Cu site), there is 
a charge gap of order e. Then the problem essentially 
reduces to that of the zigzag Hcisenbcrg chain, described 
by the effective Hamiltonian 

H e ff = J\ Sd,i,gSd,i+l,a 

— J2 2_j (Sd,i,rSd,i,2 + Sd^i-S^i+i^) (7) 
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FIG. 2: Calculated gaps from DMRG, measured in units of 
t p d, as a function of 1/L, where L is the number of Cu sites 
in the system. A c is the charge gap, A s is the spin gap, and 
A$ N =2 and As N= 4 are the gaps for transferring one and two 
holes from one chain to the other (see text). The dashed lines 
are linear extrapolations. 



Note that Ji , J 2 are different from Ji , J 2 in ([6]) — they 
arise from projecting out the p x orbitals. However, we 
still find that Ji > 0, J 2 > 0. The model (J7J) was studied 
in Ref. [lji[ by a combination of DMRG and bosoniza- 
tion, for J 2 both positive and negative. It was found 
that for J 2 > 0, the coupling between the two chains is 
irrelevant. This is due to the fact that the interaction 
between the two chains in Eq. ([7]) is geometrically frus- 
trated. Later studies found that the interaction actually 
contains a marginally relevant operato r 20 ! 21 , but the RG 
flow is extremely slow, so the system can still be regarded 
as gapless for all practical purpose 2 ^. This is an impor- 
tant difference between the zigzag ladder and the simple 
Hubbard ladder; the latter has a large spin gap A s ~ J/2 
at half filling. 



IV. DMRG SIMULATIONS 



The doped zigzag chain cannot be understood as sim- 
ply as the half filled one. Moreover, it is not entirely clear 
what doping level leads to superconductivity in Pr-247; 
counting charges (and assuming that the valence of the 
Pr ions is +3), there are 9/7 — 25/7 doped holes per cop- 
per site on average, and superconductivity occurs when 
5 > 0.3 (see Q). However, these holes are distributed 
in an unknown manner between the plane, single chain, 
and double chain Cu sites. Therefore, we have performed 
DMRG simulations with a finite doping concentration 
which we take, somewhat arbitrarily, to be 0.25 "doped 
holes" per Cu (i.e. the density of holes is taken to be 
n = 1.25 per Cu site). The following parameters were 



used in Eq. (g]): t pd = 1, e = 0.5, and J\ = J 2 = 0.5 
The values for t p d, Ji, and e arc comparable to values 
that are commonly used in effective t — J models for the 
copper-oxide planes in the cuprates. In the simulations 
described in this section, the direct oxygen-oxygen hop- 
ping, t PtP , was neglected. (The effect of non-zero t PtP is 
examined in the next subsection.) Up to m = 2700 states 
per block where kept for the longest systems (L = 88 Cu 
sites), resulting in an average truncation error of about 
10~ 6 . The energies where extrapolated to zero truncation 
errors using the standard procedure 2 ^. The convergence 
of the ground state energy as a function of the number 
of sweeps improves dramatically when a small perpen- 
dicular copper - oxygen term is added to the hopping 
Hamiltonian (J5J), with t±_ = 0.01. We have checked that 
the physical results arc not affected by adding this term. 

In Fig. [5] we plot various gaps in the spectrum as 
a function of 1/L. A c = £jv +2 + -EW-2 — 2L^v is the 
charge gap, and A s = Es= =2 — Egz=i is the spin gap. 
We define the spin gap relative to the S z = 1 state, since 
the ground state is found to have spin one for the sys- 
tems considered here^L. In addition, we have calculated 
the gap to the state with one or two holes transferred 
from one chain to the other: A$n=2 = Esn=2 — Esn=o 
and A,5jv = 4 = Egjy = 4 — Esn=o, where 5N — N 2 — N-y 
is the difference between the number of holes in the two 
chains. The direct measurement of A$n is possible due 
to the fact that Ni and N 2 (the numbers of holes on 
each chain) are separately conserved. In the presence of 
a small inter-chain hopping term, this conservation law is 
weakly violated. However, it is still possible to measure 
A^jv by applying an inter-chain potential difference AV . 
When AV = Asn, one hole is transferred from one chain 
to the other. This transition is smoothed by the inter- 
chain hopping term, with a width proportional to t± . For 
t± = 0.01, we found that this does not considerably limit 
the accuracy in the determination of Asn- 

The gaps were extrapolated linearly with 1/L to the 
thermodynamic limit (L — > 00). The charge gap A c 
extrapolates to very close to zero. The spin gap A s 
also extrapolates to a very small value, which is indis- 
tinguishable from zero to the accuracy of our calcula- 
tions. However, A,5jv=2 and Asn=a clearly extrapolate 
to finite values. The linear extrapolation results arc 
A SN = 2 (1/L -> 0) w 0.046, A 5JV=4 (1/L -»• 0) w 0.097. 

Let us assume that, as the simulation suggests, A c = 
A s = 0. In terms of the classification by the number of 
gapless bosonic modes, this implies that we must have at 
least one gapless charge mode and one gapless spin mode. 
However, since A^jy 7^ at least some of the relative 
spin and charge modes of the two chains are gapped by 
the inter-chain coupling. The fact that Asn = 4 ^ in 
the thermodynamic limit implies that the relative charge 
mode is gapped. This is because it involves transferring 
two holes from one chain to the other, and two holes can 
be combined to a singlet and therefore carry no spin. 

Fig. [3^ shows the average hole density on the oxy- 
gen sites in a chain of L = 88 Cu sites. The density 
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FIG. 3: (a) The average hole density on oxygen sites as a function of position. Blue Squares (red circles) mark the upper 
(lower) chain from the DMRG simulation with (n) = 1.25 holes per Cu site, (b) (S'f ) as a function of position. Copper sites 
are marked by squares, and oxygen sites - by circles. The size of the circles represents the average hole density on the oxygen 
sites (bigger circles mean larger hole density). 



profile shows pronounced density oscillations with a pe- 
riod of A c = 4a where a is the lattice constant. For the 
density of holes in the simulation, this corresponds to a 
wavevector of 4kp. Fourier transforming the density pro- 
file reveals that the 2kp component is almost completely 
absent, which is characteristic of systems with strong re- 
pulsion. Note also that the relative charge densities in the 
two chains appear rigidly locked to each other so that the 
charge oscillations are strictly staggered, further evidence 
of the existence of a relative charge gap. 

A zigzag chain model similar to Eq. @ was studied 
in Ref. [17| . and a spin gap was found. However, in that 
study there where no oxygen sites and it was assumed 
that J2 < (i.e. antifcrromagnctic) and |J 2 | > J\. The 
ferro mag netic case was studied (in a different context) in 
Ref. [HI j and their results seem consistent with ours (at 
least over a certain range of parameters). 

In Fig. [3Jd we show (S") along with the oxygen site 
density (represented by the size of the circles) in the same 
system. A weak Zeeman field H = O.Olz was applied to 
the first and last Cu sites on the chain in order to select a 
direction in spin space. The main periodicity in the spin 
density is X s = 8a, or 2kp. Note that a peak in the hole 
number density is always accompanied by a 7r phase shift 
in the spin density wave. Some features of the pattern 
of the charge and spin density can be understood from a 
strong coupling approach, as discussed in section V A. 



B. 



expected to be robust over a finite range of t pp . How- 
ever, this argument cannot tell us what is the range of 
stability. In order to study the effect of a finite t pp , we 
have performed simulations including an oxygen-oxygen 
hopping term H pp . Using the same notation as in Eq. 



is given by 



H n 



-t pp P 



(V ■ 



lPx,i+l,<T,2 + P X: iM,2Px,i,cr,l + h.C. ) P 

(8) 



In the presence of this term, the method we used pre- 
viously to determine whether there is a relative charge 
gap (based on measuring the energy gap for transferring 
one hole from one chain to the other) is no longer ap- 
plicable because now the charges on the two chains are 
not conserved separately. Consequently, the charge dif- 
ference cannot be controlled in the calculation. Instead, 
we used an alternative method to determine in which 
phase the system is. The local spin and charge densities 
on each chain where measured and Fourier transformed. 
We define 




(9) 



So far, we have neglected the oxygen-oxygen hopping 
t pp . This term is expected to be significantly smaller 
than tpd, because of its longer range. Estimates from 
LDA calculations^ give t pp ~ 0.27£ p <A 

As we saw, the DMRG simulation with t pp = re- 
vealed a finite gap to the relative charge mode between 
the two chains. In such a phase, single particle inter- 
chain hopping is irrelevant^. Therefore, this phase is 



Here 



n 



d,i,l 



i ad] a ^i.a^i and L is the length of the chain 
(number of Cu sites). Due to the open boundary 
conditions, n q and show pronounced peaks at cer- 
tain wavevectors. (These are known as "Friedel-like" 
oscillations^.) These peaks occur at the wavevectors 
of the gapless spin/charge modes of the system. A ID 
version of Luttinger's theorem^ states that there must 



6 



be a gaplcss mode at a certain wavevector, corresponding 
in our case to 4fcp of a single chain, where 2kp = im (n 
is the number of holes per Cu site). In addition, there 
could be other gapless modes at other wavevectors (e.g., 
charge or spin modes at 2kp for each chain, etc.). If 
there is only a single gapless charge mode, as we found 
in the case t pp = 0, then we expect to find gapless 
modes (and therefore peaks in the Fourier transformed 
local spin/charge densities of a finite system with open 
boundary conditions) only at a single wavevector, plus 
its harmonics. If, on the other hand, the relative charge 
gap closes and there is more than one gapless charge 
mode, there is no reason why these modes cannot have 
different wavevectors (these are the analogues of the 
two Fermi wavevectors in the non-interacting ladder 
problem). This way, the phase transition from a phase 
with a single gapless charge mode to a phase with two 
gapless modes can be identified. 

The results of the DMRG calculations are shown in 
Fig. [4j No additional peaks appear at new wavevectors 
and the spin/charge profiles do not change qualitatively 
until tpp rsj 0.5 — 0.6i p( j. This shows that the phase with a 
single gaplcss charge mode is robust at least up to t pp = 
0.5t pc i, which is higher than the value estimated from 
band structure calculations for the zigzag CuO chain. 
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FIG. 4: (Color online.) Fourier transforms of the charge (a) 
and spin (b) profiles (defined in Eq. ([9])) calculated by DMRG 
in a system of length L — 72 Cu sites. Various values of t pp 
(Eq. pi) were used. (t pp is measured in units of t p d-) Both 
profiles change qualitatively and new peaks appear between 
t pp = 0.5t p d and 0.6t p d, which signals a phase transition to a 
phase with two independent gapless charge modes. 



V. ANALYTICAL RESULTS IN SPECIAL 
LIMITS 

A. tpd — ^ limit 

We start from a qualitative strong coupling descrip- 
tion, in which we assume that t^ d /e -C Ji,J2- As we 
saw, at half filling the inter-chain coupling Ji is frus- 
trated, so the low-energy spectrum is essentially like that 
of two decoupled chains. A typical snapshot of the spin 
configuration in this state is depicted in Fig. [5^,. 

We now consider a lightly hole doped system. Since 
all the copper sites are occupied, the doped holes reside 
on oxygen sites. Fig. [5)3 shows two neighbor holes on 
opposite chains. Since J\ > 0, every doped hole induces 
a 7r shift in the phase of the spin density wave around 
it. The doped hole locally relieves the frustration of the 
inter-chain ferromagnetic coupling J2 , causing a net cou- 
pling between the spin fluctuations on the two chains. 
This coupling is proportional to the hole concentration 
x times J2. The hole can move along the chain via a 
second order process in t p d without disturbing the local 
spin order. Note that some exchange energy is gained 
also near the neighbor hole in the opposite chain. 

However, if two neighbor holes arc in the same chain 
(as shown in O0, the situation is different. Now we can- 
not satisfy the J 2 terms near both holes, and the inter- 
chain coupling is partially frustrated. Therefore, the low- 
est energy states will be ones in which the holes appear 
in alternating order in the two chains. The holes arc free 
to move so as to reduce their zero point kinetic energy, 



as long as they do not pass each other. 

The qualitative picture in the small t pc [ limit can ex- 
plain some of the features in Fig. [3Jd. As mentioned 
previously, the peaks in the charge density on the O sites 
are accompanied by n phase shifts in the spin density 
wave on the Cu sites, which can be understood as com- 
ing from the antifcrromagnetic coupling between neigh- 
boring Cu and O spins. The fact that the hole density 
oscillations are it phase shifted between the two chains, is 
reminiscent of the alternating order between the chains 
described above. Moreover, the spin on the O site with 
maximum local density is always parallel to the spin of 
the nearest Cu on the opposite chain, due to the ferro- 
magnetic coupling between them. The DMRG picture is 
therefore qualitatively very similar to Fig. [5] even though 
t p d is not small. 

In fact, we suspect that the present analysis is qual- 
itatively correct for small enough x independent of the 
magnitude of tpd- For a single hole in an antifcrromag- 
netic chain, a tt phase shift in the antifcrromagnetic cor- 
relations appears to be generic. Beyond that, the energy 
which favors alternating holes on the two chains is purely 
geometric in origin, and so only requires that the holes 
be sufficiently dilute. 



B. RG and Bosonization treatment 

Another limit of the zigzag chain model ([J} that can 
be treated analytically is the weak inter-chain coupling 
limit (|^2 1 Jiit)- I n this limit, only the low-energy 
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degrees of freedom of the decoupled chains are affected 
by the inter-chain coupling. These degrees of freedom 
can be described as two Luttingcr liquids. The inter- 
chain coupling can then be treated using perturbative 
RG. The nature of the strong-coupling fixed point can 
be understood using abelian bosonization. This is a stan- 
dard procedur e) 13 ' 14 . 

As we will see, in the weak inter-chain coupling limit a 
spin gap is predicted, in contradiction to DMRG results 
of section llV Al The discrepancy can be explained by the 
fact that the inter-chain coupling is not small. Rather, 
we will use bosonization to parameterize an effective low 
energy model which is consistent with the numerical re- 
sults. This model is then used to calculate low-energy 
susceptibilities of the system. 

The Hamiltonian is written as 



H = H 



LL 



(10) 



Here Hll is a Luttingcr liquid Hamiltonian that de- 
scribes the low energy degrees of freedom of each chain, 
and H int is the inter-chain coupling Hamiltonian, 



Hll — Hq s + H 



(11) 



where -ffo.s and i?o,c are bosonic Hamiltonians for the 
spin and charge modes of each chain: 

^ ^ ^ ^ ^ 

^ ^ ^ ^ ^ 

(a) 



O 



4 o \p o 4 ° t 'v4/' t 



o 



(b) 



f | f | f | $ 



(c) 



FIG. 5: Typical pattern of the spin correlations in various 
situations. Copper (oxygen) sites are shown as square (cir- 
cles), (a) The undoped chain. The intra-chain interaction is 
antiferromagnetic, while the ferromagnetic inter-chain inter- 
action is frustrated, (b) Two doped holes, one in each chain. 
Note that the phase of the spin fluctuation shifts by n around 
the doped holes. The spin of each doped holes is parallel 
to the spin of the nearest copper site in the opposite chain, 
which gains some ferromagnetic interchain exchange energy, 
(c) Two doped holes in the same chain. This is similar to (b), 
except that now the inter-chain ferromagnetic exchange en- 
ergy of the hole on the right is lost. Putting two neighboring 
doped holes in the same chain is less favorable energetically. 



H o,s = -f dx 

A 1 O " 



i=l,2 1 



K s (d x 6 hS ) 2 + {d x (p i s ) 2 



(12) 



H o,c = y Yl / dx 



i=l,2 ' 



Kc (dxBi.cf + -T7- (dxipt.cf 



and i = 1,2 is the chain index. SU(2) symmetry gives 
the constraint K s = 1, while K c < 1 depends on the 
details of the intra-chain microscopic interactions. Here 
<p and 9 represent dual fields, i.e. [d x 6 a {x), 4>b(x')\ = 
i^a,b5{x — x'), and are related to the density and current 
operators, e.g. pi = yj2/itd x 4>i, c and ji = u c y/2/nd x 8i_ c 
where pi and ji are, respectively, the electron density 
and current density on chain i. At half filling, the charge 
sector will also have a non-linear term g u cos(v / 87r</j Cj i), 
which corresponds to Akp umklapp scattering. The spin 
sector should include marginally irrelevant operators, 
which we have not written here. 

Next, we consider the interaction term between the 
two chains. This term is most conveniently written in 
terms of the corresponding fermionic degree of freedom, 
ipi,a, with i = 1,2 for the two chains, which is in turn 
written in terms of the right and left moving fermions, 
i>i,a ~ e lkFX i>iM,<j + e~ lkFX i>iX,<j, where 2k F = irn is the 
Fermi wavevector (n is the total density of holes). Taking 
a naive continuum limit of the J 2 term in ©, we get: 



H,, 



J 2 a J dx{ Si (x) ■ S 2 {x) + S 



- [ni (x) n 2 (x) + ni (x) n 2 (x + a)} } 



S 2 (x + a) 
(13) 



Here Si (x) = ip[ (x) aipi (x), rii (x) = tpj (x) tpi (x), and 
a is the lattice constant. This form satisfies the require- 
ment that the Hamiltonian is symmetric with respect to 
Si (x) — > S2 (x), Si (x) S2 (x + a) and similarly for 
the ?Vs, which corresponds to a reflection that inter- 
changes the two chains, followed by a translation by one 
Cu-0 distance. 

In bosonized form, the most relevant part of Q13p is 
written as: 



Hint — 



sin(f) 

(TTfl) 2 

91&.1 



J dxsin ^V^7riy9_ iC —J 



(14) 



Here O s \ = [cos (v^V-.s) + cos (V4tt^-,s)] and 
O s2 = [cos (%/47T( / 9_ iS ) + cos (^/iTTip +:S )] where ip± jC = 
(<£>i, c ± <P2, C ) /V2 are the even/odd charge modes of the 
two chains, and similarly for ip±_ s and 9± >s . We have de- 
fined 5 = n(n— ~ ) , where n is the density of holes per unit 
length. (<5 is a wavevector whose length is proportional 
to the amount of doping away from half filling.) At half 
filling, additional umklapp terms appear. g\ and g 2 are 
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dimcnsionlcss coupling constants, whose bare values at 
the initial scale are: g\ = —g^ = Under RG, g\ and 
172 will flow, and need not remain equal in magnitude. 
However, additional couplings are prevented as long as 
the exact SU(2) spin-rotational symmetry is respected. 
Specifically, there are three distinct cosines of spin fields, 
but only two independent coupling constants. 

Note that the inter-chain coupling (fT3|) vanishes in the 
limit 5 — > 0, i.e. at half filling. This is due to the frustra- 
tion of the inter-chain coupling in that limit. (See section 
IIIII ) Then marginal operators need to be considered. 
These produce an extremely small gap in the spectrum, 
especially in the J2 < case^, so the system can be 
considered as essentially gapless at half fillin g 19 i 20 . 

Away from half filling, both the g\ and 172 terms in (fT4|) 
are relevant, and produce a gap in the spectrum. They 
both have the same scaling dimension of 1 + K+ )C < 2 
(assuming that K± iS = 1). The system flows to a strong 
coupling fixed point at which the <^- )C , <p+, s and 
fields are pinned. Then the only gapless mode is the 
total charge mode, and the system is in a C1S0 phase, 
similar to the generic situation in the Hubbard ladder—. 

However, the DMRG result indicates that while there 
is a substantial gap in the relative charge mode, the gap 
in the total spin sector is either zero or very small (see 
Fig. [2]). This discrepancy is likely to be caused by less rel- 
evant (sub- leading) operators, not included in Eq. (TH)) . 
whose bare coefficients are of order unity (since we are 
not initially at weak coupling). The neglect of these op- 
erators in the initial stages of the RG transformation is 
a quantitatively unreliable approximation. Taking the 
DMRG result into account, we hypothesize that the ef- 
fective value of gi at the fixed point is close to (since 
non-zero <?2 is what induces a full spin-gap). If this is 
true, then over a broad intermediate range of energies, 
the system is governed by the unstable fixed point with 
92 = and non-zero g\. 

We would like to stress that the smallness of the spin 
gap is not a result of fine tuning, but rather appears 
(from the numerics) to be a ubiquitous property of the 
model. (The same result is obtained over a range of dif- 
ferent parameters and doping levels.) In contrast, if the 
inter-chain coupling J2 is artificially turned to be antifer- 
romagnctic, a spin gap appears in the numerical simula- 
tion. Therefore the small spin gap seems to be a feature 
of the ferromagnetic J2 case. 

Note that the phase with 171 ^ 0, gi = is unusual, 
because neither the field (f- >s nor its dual #_. s can be 
pinned (since O s \ contains the cosine terms containing 
both fields with equal magnitude). The properties of this 
phase, which follow closely from an earlier analysis of the 
Heisenbcrg ladder by Shclton, Nersessyan and Tsvelik^S, 
are explored in appendix [A] The main conclusion is that 
essentially "half of the relative spin mode is gapped, and 
the other half is gapless. Therefore we denote this phase 
as C1S3 . 



C. Physical susceptibilities 

Given that the system is in a C1S3 phase, the physical 
susceptibilities can be calculated in a similar way to that 
described in [26[, with addition of the charge modes. For 
an example of such a calculation, see appendix El The 
leading temperature dependence of these susceptibilities 
depends only on the total charge Luttinger parameter 
A' + , c . The superconducting susceptibility is 

Xsc-T^fe" 1 (15) 

for both singlet and triplet pairing, which gives that the 
superconducting susceptibility is divergent for K + c > 
0.4 (i.e. even for strongly repulsive interactions). 
Other susceptibilities are 2&f SDW and CDW, with 

K +,c 5 

XSDW,2k F ~ XCDW,2k F ~ T 2 * (16) 

and 4fc F CDW: 

XCDWAk F ~ T 2K+ ' C ~' 2 (17) 

The superconducting susceptibility is the most divergent 
one for K+ )C > 1. 

Assuming that the relative charge mode tp— tC IS 
gapped, it is possible to extract K + c from the den- 
sity profile in the DMRG simulation (sec Fig. [3]). 
The density profile is expected to behave as (rij) ~ 
cos(AkpXi)/[Lsm(TTXi/L)] K +-°, which is essentially the 
square root of the Akp part of the density-density corre- 
lation functional. The amplitude of the Akp CDW near 
the middle of the chain, A Xi —L/2 : thus decays as L~ K +>°. 
Fitting A Xi=L j 2 of chains with different lengths to this 
expression, we obtain K + . c = 0.6 ± 0.05. 

Since NQR measures of the spin-lattice relaxation rate 
have been done on superconducting Pr-247, it is interest- 
ing to extract the temperature dependence of this quan- 
tity from the theory. Assuming that the main relaxation 
mechanism is the coupling of the nuclear spins to fluctu- 
ations of the local magnetic field due to electronic spins, 
1/Ti is given by-21: 

Y 1 =^J2\ A ^i)\ 2 s aa ( q ,i, R ) (18) 

qa 

where A a (q) is the form factor of the hyperfine Hamil- 
tonian, cor is the nuclear spin resonance frequency, and 
S aa (q, LOji) is the electronic spin structure factor, ljr 
is typically smaller than any other energy scale in the 
problem, so we will assume l>r = 0. 

The spin structure factor in the CiS| is expected to 
be of the form: 

S aa (q,Lu = 0)~ Alnq + Bq^-i (19) 

where A and B are constants. The two terms here are 
the two main contributions to S aa (q,u> = 0) which come 
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the vicinity of the points q = 0,2kp. Integrating (fH 
over q, we get the dominant temperature dependence of 
the spin-lattice relaxation rate: 



1 

ATlnT + BT— 



(20) 



NQR measurements^ show that T^ 1 behaves as a power 
law of temperature, with different exponents above and 
below T c . The fact that a power law is seen above T c 
is a clear evidence for the absence of a spin gap, and is 
consistent with a C1S3 phase. According to the mea- 
surement, -t- oc T 5 above T c . Comparing this behavior 
to (|20|) . we see that it corresponds to K+. c — 1.5. This is 
consistent with dominant superconducting correlations. 

The value of K + c extracted from the NQR data is sub- 
stantially larger than the one extracted from our DMRG 
calculation. Clearly, this is a significant discrepancy. 
However, it is well known that K +tC is a non-universal 
exponent, and in this case it depends strongly on micro- 
scopic details of the problem. Since the detailed aspects 
of the microscopic model are certainly not "realistic," we 
do not feel that any microscopic calculation can be ex- 
pected to predict the experimental value of K + c reliably. 
A better strategy is thus to extract K + c directly from 
experiments. As we saw, in the superconducting sample 
this yields a value of K +iC that corresponds to dominant 
superconducting fluctuations. 

The NQR measurement was done on with a sample 
with an oxygen content S = 0.5 which is close to the 
optimal value for superconductivity. We are not aware 
of any similar NQR data on samples with a different O 
content and T c . Were such data obtained, we would pre- 
dict a smaller power should govern the T dependence of 
l/Ti(T), corresponding to a lower K +)C . 



VI. CONCLUSIONS 

To summarize, we have considered the possi- 
bility of zigzag chain-driven superconductivity in 
Pr2Ba4Cu70i5_5. Assuming that the chains are weakly 
coupled, this implies that the single zigzag chain must 
have a large superconducting susceptibility. This, in con- 
junction with the fact that Cu NQR experiments do not 
show any spin gap above T c , raises the possibility that 
the zigzag chain is in a phase in which some, but not 
all, of the modes of the two-component electron gas are 
gapped. We have shown evidence for such a phase using 
a microscopic copper-oxygen model for the zigzag chain. 

The gapping of some of the relative spin and charge 
modes enhances superconducting (as well as other) fluc- 
tuations at low temperatures. The pairing operator is 
composed of one hole from each chain with zero center 
of mass momentum. Since there is no spin gap, triplet 
and singlet superconducting fluctuations are equally en- 
hanced. 
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APPENDIX A: THE C1S3 PHASE 

2 

We would like to describe a phase of the two- 
component electron gas with a gap in the relative spin 
and charge sectors, but no gap in the total spin and 
charge sectors. This is the phase that seems to emerge 
from the DMRG results. (See Fig. H) Starting from 
the Hamiltonian H = Hll + -ffmt, where Hll is a 
Luttinger liquid intra-chain Hamiltonian pip and H int 
is the most relevant part of a generic SU(2) invariant 
inter-chain interaction (Eq. 1 14|) . we see that in or- 
der for the spin gap to vanish we must have 52 = 
in (|14p ). The remaining term is then proportional to 
O s \ = cos (v47n/3-, s ) +cos (\/47r0-,s). The total scaling 
dimension of this term with respect to the Luttinger liq- 
uid fixed point is 1 — A'_ iC , with A'_. c < 1, so it will grow 
under an RG transformation. At some scale, g\ becomes 
of the order of unity, and we can replace sin (\/47r<yS_. c ) 
and cos (\/47r(yS_. c ) by their mean value (since the field 
<y2_. c will be pinned to the minimum of the potential). 
However, since the interaction contains the cosines of the 
conjugate fields </?_, s and 9- tS with an equal weights, both 
fields fluctuate strongly at the fixed point, and neither is 
pinned. Wc arc faced with the problem of how to describe 
the low-energy properties of the resulting phase. 

A remarkably elegant solution for this problem is given 
in a paper by Shclton et al [2(j. They considered the 
problem of two weakly coupled spin 1 chains, but the 
results are readily generalized to our case by neglecting 
the fluctuations in the c field, replacing it by its mean 
value. The Hamiltonian (TTUJ) can then be solved by re- 
fcrmionizing of the fields tp± tS - We introduce two Dirac 
fermions, ifi and x, as follows: 

V27ra 



1 



XR,L = 



tv^f(9-,»±(p-,,) 



(Al) 



where R and L denote right or left moving fields. The 
Hamiltonian for -0 is then just a free Dirac Hamiltonian: 



H + = u 



The interaction Hamiltonian Hi, 



(A2) 



becomes 
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quadratic in x, giving 

H - = Us J dx { x «ii XR ~ x ^ XL ) 

+ ^ Jdx (xW + xWr H.c.) (A3) 

Here M depends on g\ and the average of the f- lC part 
in (|13[) . (|A3|) is most conveniently diagonalizcd by writ- 
ing x m terms of two Majorana fermions. Adopting the 
notation of [291 , we write 



Pu 



~V2~ 



Xv - XI 
iV2 



(A4) 



with v = R, L. Plugging this into (|A3[) . we get: 

H- = H [f] + F M W (A5) 



Id Id" 

-~tpr - pl- -j- pl 

i dx i ax 



■ [ dx (pr--^- 

J \ i ax 

l — J dx {p R pL - PlPr) 



(A7) 



The field £ is thus massless, while p is massive, with 
a mass M. In that sense, "half of the field f-, s is 
gapped. The total spin sector, on the other hand, is 
completely massless. We therefore denote this phase as 
a C1S3 phase. 

Setting 52 7^ in (fT3"|) will gap out also the total spin 
mode, giving a C1S0 phase. Moreover, adding less rel- 
evant inter- or intra-chain operators to (|13p would gen- 
erate the <72 term, even if it is not present in the bare 
Hamiltonian. However, the DMRG calculation presented 
in section II V A I indicates that the spin gap A s is much 
smaller than Asn, the gap to transferring one hole from 
one chain to the other, which is related to a gap in 
the relative spin/charge modes. Therefore, at interme- 
diate energies (or temperatures) between A s and A^jy, 
the physics is expected to be dominated by the unstable 
C1S3 fixed point. 

Next, let us find the long-distance behavior of the 
correlation functions of physical operators in the C1S3 
phase. For concreteness, we will focus on the pair field 
operator Asc = ^ia.t^l.j.- This operator creates a 
pair of holes in the two chains with opposite momenta 
and spins. In bosonized form, this operator becomes 



9 +c is a free bosonic field, whose long-range correla- 
tions are (e i ^ e +^e l ^ +-^\ 
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Similarly, 



x 2 5 since (p+, s is a free field 

with K+,s — 1, as dictated by the SU(2) symmetry. This 
correlation function is expected to have logarithmic cor- 
rections due to marginal operators, which we have ne- 
glected. The <p~ c held is massive, so at long distances 
we can replace e'-V^v-.c by its expectation value. The 
treatment of e l ^ e_,s is more subtle, since as we have 
mentioned, this field is "half gapped" in the C1S3 phase. 
It is nevertheless possible to calculate its correlations, fol- 
lowing a method used in [2(|. We describe this method 
briefly here. The relative spin sector, which is described 
as two independent Majorana theories, can be further 
mapped onto two 2d Ising models. Since one of the Ma- 
jorana fields is massless and the other is massive, one 
of the corresponding Ising fields is at criticality and the 
other is away from criticality. It can be shown that the 
four operators sin (^/n6- tS ), cos (t/tt6- )S ), sin (^/n(p- tS ), 
cos (y / 7T(/5_ s ) have the following corrcspondance to the 
order and disorder operators of the two Ising modcl o 26 ' 28 : 



cos s/n(p- ia 
cos ^/ttO-^s 



P1P2, 



sin v / 7T(/3_ i 
sin ypn®- i 



- (710-2 

01M2 (A9) 



Here 01, p,\, cr 2 , p.2 are the order/disorder operators of 
the two Ising models labelled 1 and 2. Note that the 
scaling dimension of all the operators in the left hand side 
of (|A9[) at criticality is |, which is consistent with the fact 
that the dimension of the Ising operators at criticality is 
g. Carrying out the correspondence to the Ising model 
carefully, one finds that the Ising model labelled by 1 
is in its disordered phase, so (p,x) ^ 0, (01) = 0. The 
other Ising model labelled by 2 is critical, and therefore 
(<72 (x) 02 (0)) ~ x~ 1 , and similarly for p.2- Thus at long 
distances, ^ e iV^9-A^) e iV^9-, 3 (o)\ _ x -{_ The exponent 

of this correlation function is just half of the exponent 
that one gets for a free 9- tS field. This can be understood 
as a consequence of the fact that "half of the — , s mode is 
gapped. Therefore, correlation function of Asc behaves 
as 



A^ (x) A sc (0) 



1 



(A10) 



A sc - e l ^ {e ^ +e - 



S +V-,c + <P+, S ) 



(A8) 



which gives a superconducting susceptibility that behaves 

1 5 

as xsc (T) ~ T 2K +> C 4 . The correlation functions of 
other physical operators (CDW, SDW etc.) can be cal- 
culated in a similar manner. 

It is interesting to note that the Hamiltonian (JX5J) with 
g-2 = is equivalent to the integrable super-symmetric 
super-Sinc-Gordon model 2 ^. 
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